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O '■ ABSTRACT 

23: 

We have used high resolution WFPC2-HST and ground based Wide Field 
images to determine the center of gravity and construct an extended radial den- 



sity and brightness profile of the cluster NGC 6752 including, for the first time, 
detailed star counts in the very inner region. The barycenter of the 9 innermost 
X-ray sources detected by Chandra is located only l'/9 off the new center of grav- 
ity. Both the density and the brightness profile of the central region are best 
fitted by a double King model, suggesting that NGC 6752 is experiencing a post- 
core collapse "bounce" . Taking advantage from our new optical data, we discuss 
the puzzling nature of the accelerations displayed by the innermost millisecond 
pulsars detected in this cluster. We discuss two possible origins to the accelera- 
tions: 1) the overall cluster gravitational potential which would require a central 
projected mass to light ratio of order 6-7 and the existence of a few thousand 
solar masses of low-luminosity matter within the inner 0.08 pc of NGC 6752; 2) 
the existence of a local perturber(s) of the pulsar dynamics, such as a recently 
proposed binary black hole of intermediate (100-200 M & ) mass. 



Subject headings: Globular clusters: individual (NGC 6752); stars: evolution - 
binaries: close; neutron stars: millisecond pulsar 
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1. Introduction 

Although NGC 6752 is one of the nearest clusters, there is no consensus in the literature 
on its dynamical status. It was initially classified as a post-core collapse (PCC) by Djorgovski 
& King (1986, hereafter DK86) and Auriere & Ortolani (1989, hereafter A089). Later 
Lugger, Cohn & Grindlay (1995, hereafter LCG95) argued that the radial profile was not 
inconsistent with a King model. As part of a project devoted to the study the characteristics 
of globulars harboring millisecond pulsars (MSPs), we present a new determination of the 
center of gravity (§3) and of the star density profile (§4) of NGC 6752 both based on new 
high-resolution and wide-field observations. By coupling the superior resolution of HST 
with the imaging capability of the wide field camera (WFI) mounted at the 2.2m ESO-MPI 
telescope, we have obtained the most accurate and extended radial profile ever published for 
this cluster. After modeling these data (§5) we examine the dynamical status of the cluster 
(§6). 

Recently, observations of \P / P\ and location of five millisecond pulsars in NGC 6752 
(D'Amico et al. 2002) have suggested a surprisingly high mass-to-light ratio Ai/Cy in its 
core and the occurrence of non thermal dynamics in the inner regions (see Colpi, Possenti & 
Gualandris 2002). The determination of M/Cy by D'Amico et al. relied on results of pulsar 
timing and published optical data derived from medium resolution ground based observation 
only. Taking advantage from our new optical data, we re-examine in §7 the D'Amico et al. 
measurement, discussing the possibile origin and the consequences of the observed values of 
\P/P\ at some length. 



2. Observations and data analysis 

The photometric data used here consist of two sets: (i) — High resolution set — a series 
of HST- WFPC2 images was obtained in March 2001, through the F555W (V) and F336W 
(U) filters as part of a long term project (GO-8709, PI: F. R. Ferraro) aimed to study the 
central stellar populations in a set of Galactic Globular dusters (GGCs). In this dataset the 
planetary camera (PC, which has the highest resolution ~ O'.'046/pixel) was roughly centered 
on the cluster center while the Wide Field (WF) cameras (at lower resolution ~ O'.'l/pixel) 
sampled the surrounding outer regions; (ii) — Wide Field set — a complementary set of multi- 



1 Based on observations with the NASA/ESA HST, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. Also based on WFI observations 
collected at the European Southern Observatory, La Silla, Chile, within the observing programmes 62. L- 
0354 and 64.L-0439. 
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filter (B, V, I) wide field images was secured during an observing run at the 2.2m ESO-MPI 
telescope at ESO (La Silla) in July 1999, using the Wide Field Imager (WFI). The WFI has 
exceptional imaging capabilities — the image consists of a mosaic of 8 CCD chips (each with 
a field of view of 8' x 16') giving a global field of view of 33' x 34'. The cluster was roughly 
centered on chip #2. 

2.1. Photometry 

The raw WFI images were corrected for bias and flat field, and the over-scan region was 
trimmed using standard IRAF 2 tools. The point spread function (PSF) fitting procedure was 
performed independently on each V and B images, using DAOPHOT II (Stetson 1994). A 
final catalog listing the instrumental B, V magnitudes for all the stars in each field has been 
obtained by cross-correlating the B and V catalogs. The WFI catalog was finally calibrated 
by using the data-set of Buonanno et al. (1986). 

The photometric reductions of the high resolution images were carried out using RO- 
MAFOT (Buonanno et al. 1983), a package developed to perform accurate photometry 
in crowded fields and specifically optimized to handle under-sampled PSF (Buonanno & 
Iannicola 1989) as in the case of the HST-WF chips. 

PSF-fitting instrumental magnitudes have been obtained using the standard procedure 
described in Ferraro et al. (1997a, 2001). The final catalog of the F555W and F336W 
magnitudes was calibrated by using the zero-points listed by Holtzman et al. (1995). 

2.2. Astrometry 

The recently released Guide Star Catalog (GSCII) was used to search for astrometric 
standards in the entire WFI image field of view. Several hundred astrometric GSCII refer- 
ence stars were found in each chip, allowing an accurate absolute positioning of the sources. 
An astrometric solution has been obtained for each of the 8 WFI chips independently, by 
using a procedure developed at the Bologna Observatory (Montegriffo et al. 2002, in prepa- 
ration). At the end of the entire procedure, the rms residuals were of the order of ~ 0'.'2 
both in RA and Dec. 



2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Fundation. 
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The small field (2(5 on the side) of the high resolution WFPC2/HST images was entirely 
contained within the field of view of the WFI chip #2, which imaged the central part of the 
cluster. Thus we used more than 1200 bright stars in the WFI catalog lying in the WFPC2- 
FoV as secondary astrometric standards in order to properly find an astrometric solution for 
the WFPC2 catalog. We estimate that the global uncertainty in the astrometric procedure is 
less than ~ 0"4 both in RA and Dec. At the end of the procedure the two catalogs (WFPC2 
and WFI) have a fully homogeneous absolute coordinate system. 

The two lists were then matched together: stars in the overlapping area were used to 
recheck the homogeneity of the calibration of the ^-magnitude (which is used to construct the 
surface profile) in the two catalogs. To do this the HST F555W band was transformed into the 
Johnson V system and a homogeneous list of V magnitudes and absolute coordinates for the 
sources in the HST and the WFI catalogs was produced. In order to avoid strong selection 
effects due to possible incompleteness of the samples (in particular the WFI catalog), we 
grouped the two samples as follows: (1 ) HST sample: all stars in the WFPC2 FoV with 
V < 18.5 and r < 96" from the cluster center; (2) WFI sample: all the stars in the WFI 
FoV with V < 18.5 and r > 120" from the cluster center. The Color-Magnitude Diagrams 
(CMDs) derived from these samples are shown in Figure 1. 



3. The center of gravity 

The first step toward the computation of the density profile is the determination of 
the center of gravity C grav of the cluster. In doing this, we estimated the position of the 
geometrical center of the star distribution taking advantage from the knowledge of the 
exact star positions even in the innermost central region. We applied the procedure de- 
scribed in Montegriffo et al. (1995) which computed C grav by simply averaging the a and 
5 coordinates of stars lying in the PC camera of the HST catalog. C grav is located at 
aj2ooo = 19 h 10 m 52:04, 5 J20 oo = -59°59'04"64 with a typical la uncertainty of 0"5 in both 
«J2ooo and 5/2ooo> corresponding to about 10 pixels in the WFPC2/HST images. 

A089 derived the C grav by using the barycenter of the bright (V < 16) resolved stars 
from ground based observations. We transformed the center position in their Figure 3, 3 into 
our coordinate system. It well agrees with our determination, being located only ~ 2" S. 

The coordinates of the cluster center previously reported in the literature refer mostly 



3 Note that the orientation in both the maps shown by A089 — Figures 1 and 3, respectively — is incorrectly 
reported in their captions since E is on the right and N is at the bottom. 
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to the center of luminosity C\ um , which is usually determined by the so-called mirror auto- 
correlation (see Djorgovski 1988). C grav turns out to be ~ 10" S and ~ 2" E of the C\ UTa 
reported in the Djorgovski (1993) compilation. Such a difference is not surprising — similar 
offsets have been found in other clusters (see for example 47 Tucanae — Calzetti et al. 1993, 
Montegriffo et al. 1995; M80— Ferraro et al. 1999a). 

Figure 2 shows the central 20" x 20" of the cluster with respect to the C grav determined 
in this work (marked in Figure with a large cross at (0,0)). The center position by A089 
(small cross labeled with AO) and by Djorgovski (1993; small cross labeled with D) are also 
indicated. 

Pooley et al. (2002) noticed that the highest concentration of X-ray sources detected by 
Chandra in a 30 ksec pointing toward NGC 6752 appeared surprisingly displaced with respect 
to the optical center by Djorgovski (1988). On the contrary, Figure 2 shows that our value of 
Cg rav lies close to the barycenter of the nine more central objects (marked in Figure 2 with a 
large cross labeled with XRS at coordinates (a J20 oo = 19 h 10 m 5r!82, <5 J20 oo = -59° 59' 03"84), 
reconciling their projected positions with an almost spherical distribution. 

4. The surface brightness and star density profiles 

By using the combined data set shown in Figure 1 we computed star density and surface 
brightness profiles applying the standard procedure fully described in Calzetti et al. 1993 
(see also Ferraro et al. 1999a, Paltrinieri et al. 2001). The entire photometric sample has 
been divided in 35 concentric annuli centered on C grav , spanning a spatial range from 0" to 
27'; each annulus has been in turn split in a number of sub-sectors (generally octants or 
quadrants depending on the shape of the FoV covered by the HST and WFI fields). The 
surface brightness of each sub-sector has been evaluated as the sum of the luminosity of all 
the stars lying in it, divided by the area (expressed in arcsec 2 ). The average and the rms 
over the brightness of the sub-sectors determine the value and the uncertainty of the surface 
brightness of an annulus. However, evaluating the brightness over small regions can suffer by 
large fluctuations due to small number statistics of few bright giants. Hence, we computed 
three radial profiles, removing the stars brighter than V = 12, V = 13, V = 14, respectively. 
Figure 3 shows that the overall structure of the profile does not change with the adopted 
magnitude limit. In contrast, the rms between sub-sectors significantly decreases once the 
brightest stars are excluded. The brightness profile by A089 is shown (as open circles) for 
comparison in the lower panel of Figure 3. The agreement with the results of this work is 
excellent, although their profile is slightly fainter in the very central region, probably due to 
the incompleteness of their sample. This might be due to the fact that removal of the giants 
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can lead to an unintentional subtraction of the contribution of close, fainter stars, which are 
often unresolved in ground-based observations. 

We can take advantage from the fact that we resolve the stars even in the innermost 
region of the cluster and directly produce a star density profile. As noted by LCG95, a 
star density profile is not affected by the statistical fluctuations caused by the bright stars, 
and it is the most robust mean for determining the cluster structure parameters. Once the 
photometric sample has been divided in annuli and sub-sectors, the number of stars lying 
within each sub-sector has been counted and averaged. Star density has been eventually 
obtained by dividing the average star number by the area of the sub-sector (expressed in 
arcsec 2 ). The star density profile values obtained from this procedure are listed in Table 1. 
This is the most complete and extended density profile ever published for this cluster, since 
it samples the cluster population from the very inner region out to r ~ 27' from the cluster 
center. 



5. Modeling the radial density profile 

Usually GGCs are considered core-collapsed or not depending on how well their radial 
distribution of stars is fitted by King (1966) profiles (Trager et al. 1995). These models are 
characterized by two parameters, the core radius, r c , and the tidal radius, r t , or, alternatively, 
the concentration c = log(r t /r c ). As stated by Meylan & Heggie (1997, hereafter MH97) 
a more general rule is that all clusters with a concentration parameter c > 2-2.5 can be 
considered as collapsed, on the verge of collapsing or just beyond the collapse. 

In order to reproduce the observed profile, we obtained the surface density by projecting 
the star density from a standard isotropic, single-mass King-model derived using the code 
described in Sigurdsson & Phinney (1995). 4 The result is shown in Figure 4. As can be 
seen, a King model which properly fits the outer region significantly overestimates the star 
density over an extended inner region (5-20") . This fit closely resembles the result obtained 
by LCG95 who fitted the [/-band surface brightness profile for r < 120" with a c = 2 King 
model (see their Figure 2r). On the basis of that fit, they concluded that the cluster is not 
required to be in a post-core collapse state because the model did not differ from the data in 
a statistically significant way. However, there were aspects of the fit which led us to explore 
other possibilities. These included the bad fit in the region 5 < r < 20" and the fact that 
even in the outer parts there are regions where the data line systematically above or below 



4 In our earlier work we have used an analytic representation of King profiles. Parameter values derived 
from the analytic fits differ somewhat from those derived using the current scheme. 
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the observed points. We suspected that unmodeled effects were degrading this fit. 

For this reason, we searched for alternative solutions. As shown in Figure 5 the density 
profile can be well fit by a dual King model. The outer cluster (r > 10") is well fit by a model 
with r c = 28" and c = 1.9. The observed counts are significantly in excess of this model 
within the central 8". In this innermost region a King model with c = 2.1 and r c = 5.7" fits 
the data, but lies significantly below all of the points with r > 10". A dual King model does 
not represent a detailed equilibrium of cluster structure. However, it is consistent with the 
scenario in which the central regions have evolved away from a global King model that now 
only characterizes the outer regions. 

This anomalous structure of the innermost radial profile of NGC 6752 has been evident 
since the very first studies (Da Costa 1979). The [/-band data presented by DK86 showed 
a sharp shoulder in the radial profile and prompted the authors to classify NGC 6752 as a 
core-collapsed cluster. A few years later, A089 modeled their detailed U-band profile with 
a power law in the region between 3" and 60" from the center, noting a significant flattening 
for r < 3". 

LCG95 used also a modified power-law (see their eq. 3, also reported in Figure 6), 
in order to test the presence of resolved cores. We used the inner part of the star density 
profile of Figure 5 to repeat an analysis similar to that presented by LCG95. Besides a 
normalization factor, their expression is a function of two free parameters: the power-law 
index a and the scale- length of the core r , which can be easily converted in the usual r c 
(see their eq. 2). As can be seen in Figure 6 a proper fit (with x 2 /d.o.f = 0.28 for 8 d.o.f) 
can be obtained, yielding tq = 3"1 ± I'M (corresponding to a core radius r c = 5"2 ± 2"4) 
and a = —1.05 ± 0.05 (lex errors). These values are consistent with the previous results by 
LCG95 (r c = 6"7 ± 1"9 and a = -0.97 ± 0.15) and by A089, who estimated a radius of 
the central plateau of ~ 3" (from the analysis of their images) and a = —0.95 (from the 
brightness density profile). At a distance of D = 4.3 ± 0.4 kpc (Ferraro et al. 1999b), the 
core radius we inferred corresponds to a physical dimension of 0.11 ± 0.05 pc. 



6. The dynamical status of NGC 6752 

The significant deviation of the star number density profile from a canonical King model 
is a clear indication that the innermost region NGC 6752 has experienced (or is experienc- 
ing) a collapse phase. Unfortunately, no similarly unambiguous signature is available for 
differentiating the in- and post-core collapse state (Meylan & Heggie 1997). However, in 
favorable cases the phase of the collapse can be evaluated from indirect evidence. 
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Our results solidify the earlier suggestions of post-core collapse bounce made by A089. 
Post-core collapse clusters are expected to undergo large amplitude oscillations in core size 
due to the gravothermal instability of collisional systems (e.g. Cohn et al. 1991). The 
oscillating core spends most of the time at near maximum size and a radial extension of 0.11 
pc is consistent with the maximum radius of the core predicted by the models of post-core 
collapse bounce. The parameter most commonly used in theoretical studies is the ratio of 
core radius to half-mass radius, r c /r h . Using our value of r c = 5.2" and the half-mass radius 
of 115" from Trager et al (1995), we find a ratio r c /r h = 0.045. For comparison, the multi- 
mass Fokker-Planck models for the post-collapse evolution of M15 presented by Grabhorn et 
al (1992, see their Fig. 5) reach a comparable value of r c /r h = 0.034 during the maximally 
expanded state of the most extreme core bounces. The inclusion of primordial binaries in 
Monte Carlo simulations by Fregeau et al (2003) results in an even wider range of predicted 
r c /r h in the post-collapse phase. 

Furthermore, we have found that there is an intermediate region in the composite density 
profile of NGC 6752 which is well represented by a power law profile with a slope a ~ — 1, 
compatible with the steepness predicted by single mass models of expanding bouncing cores 
(LCG95). 

It is also worth noting that NGC 6752 has apparently retained a substantial primordial 
binary population (Rubenstein & Bailyn, 1997); these binaries may play an important role 
in supporting the core and delaying the core-collapse event. In this respect, Ferraro et al 
(1999a) have suggested that some species originated from binary evolution could be used as 
possible tracers of the cluster dynamical evolution. In particular, the large BSS population 
recently found in M80 by Ferraro et al (1999a) might be the signature of a transient dynamical 
state, during which stellar interactions are delaying the core-collapse process leading to an 
exceptionally large population of collisional-BSS. On the other hand the BSS population 
found in the central region of NGC 6752 is small (Ferraro et al 2003, in preparation) perhaps 
indicating that NGC 6752 is in a different dynamical evolutionary state than M80. Maybe 
the binary population in NGC 6752 has not been "burned out" producing collisional BSS 
while that in M80 has. 



7. The interpretation of the MSPs accelerations 

NGC 6752 hosts 5 known millisecond pulsars (D'Amico et al. 2001a; D'Amico et al. 
2002). The positions in the plane of the sky of three of them (PSRs B, D and E, all isolated 
pulsars) are close to the cluster center, as expected on the basis of mass segregation in 
the cluster. In particular PSR-B and E show large negative values of P, implying that the 
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pulsars are experiencing an acceleration with a line-of-sight component a/ directed toward the 
observer and a magnitude significantly larger than the positive component of P due to the 
intrinsic pulsar spin-down (see e.g. Phinney 1992). What is the origin of such acceleration? 
Given the location of NGC 6752 in the galactic halo and knowing its proper motion (Dinescu 
et al. 1999), it is possible to calculate the contributions to P due to centrifugal acceleration, 
differential galactic rotation and vertical acceleration in the Galactic potential, all of them 
resulting negligible (D'Amico et al. 2002). Hence, the remaining plausible explanations of 
the observed negative P are: the accelerating effect of the cluster gravitational potential well 
or the presence of some close perturbator(s) exerting a gravitational pull onto the pulsars. 
Taking advantage from the new results presented in this paper we discuss in the follow 
viability and implications of these two possibilities. 

7.1. Case (i): Overall effect of the GC potential well 

The hypothesis that the line-of-sight acceleration of the MSPs with negative P is dom- 
inated by the cluster gravitational potential has been routinely applied to many globulars. 
In particular from this assumption a lower limit to the mean projected mass-to-light ratio 
in the U-band Ai/Cy in the central region of M15 (Phinney 1992) and 47 Tucanae (Freire 
et al. 2003), yielded Ai/£y> 3 and > 0.7 respectively. Following Phinney (1992), a lower 
limit to Ai/Cy in the inner regions of NGC 6752 is given by 

ley 

- 1-1? = 5-1 x 10-£ (J^) s"\ 

where Sy(< 9±) is the mean surface brightness within a line of sight subtended by an angle 
9j_ with respect to the cluster center and M cy i{< 9±) is the mass enclosed in the cylindrical 
volume of radius R± = D9j_. This equation holds to within ~10% in all plausible cluster 
models and is independent of cluster distance, except for the effects of extinction. Since 
E(B - V) is very small for NGC 6752 (= 0.04 according to Harris 1996) the latter is a 
negligibe affect for this cluster. 

Using the observed P/P of PSR-B and PSR-E (D'Amico et al. 2002), combined with 
our accurate determinations of C grav (the cluster center of gravity) and £y(< r) (the mean 
surface brightness radial profile), it turns out that Ai/Cy > 6-7 (see Figure 7) for the case 
of NGC 6752. D'Amico et al. (2002) obtained a slightly larger Ai/Cy > 10 using published 
values of C grav and Sy (< r) derived from medium resolution ground based observations only. 



< 
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The difference between the two estimates is mainly due to our new position of the cluster 
center of gravity. Despite a residual uncertainty ~ Of! 7 on C grav , under the hypothesis that the 
line-of-sight acceleration of PSR-B and PSR-E are entirely due to the cluster gravitational 
potential, a lower limit of AXjLy > 5.5 can be firmly established. It is obtained assuming 
that the two millisecond pulsars were just symmetrically located (and hence at the minimum 
projected distance) with respect to the actual center of gravity. 

The sample of the core collapsed clusters shows typical values of the projected central 
mass-to-light ratio in the interval 2-3.5 (Pryor & Meylan 1993), although larger Ai/C ratios 
can be obtained when a Fokker-Planck model fit is used (e.g. the case of M15 Dull et al. 1997, 
2003). If we take Ai/Cy ~ 3 the expected total mass located within the inner t^b = 0.08 
pc of NGC 6752 (equivalent to the projected displacement of PSR-B from C grav ) would be 
~ 1200-2000 M Q . On the other hand, the observed Ai/Cy^ 6-7 implies the existence of 
further ~ 1500-2000 M of low-luminosity matter segregated inside the projected radius 
r This extra amount of mass could be constituted by a relatively massive black hole, like 
the ~ 17001?™ M & black hole in the center of the globular cluster M15 recently proposed 
by Gerssen et al. (2003). However, the results presented in §4 and §5 show that there are 
significant differences in the kinematics and mass distributions of the central regions of M15 
and NGC 6752 (we note here that because the distance to NGC 6752 is less than half that 
to M15, its inner core is more easily studied): (1.) HST imaging of M15 (Guhathakurta 
et al. 1996 — WFPC2; Sosin & King 1997 — FOC) shows no evidence for a compact core, at 
variance with our observations of the core of NGC 6752; (2.) The derived stellar density 
profiles of M15 have power law slopes consistent with a = —0.75 (expected from single mass 
models with a dominant central black hole; LCG95). If a ^ 1O 3 M black hole resides in the 
inner region of NGC 6752, its gravitational influence would extend more than ~ 2" from 
the center of the cluster and probably produce a central power-law cusp, which we do not 
observe. 

A very high Ai/£y ~ 6-7 could be also due to central concentration of dark remnants 
of stellar evolution, like neutron stars (NSs) and heavy ~ 1.0 M white dwarfs (WDs) (as 
also proposed for M15 by Dull et al. 1997, 2003 and by Baumgardt et al. 2003). In this case, 
one can constrain the initial mass function (IMF) and/or the neutron star retention fraction 
/ rct in NGC 6752. Based on the current population of turn-off stars (in the mass interval 
0.6-0.8 Mq), the estimated number of upper main sequence initially present in NGC 6752 
is ~ 4000 (D'Amico et al. 2002). This assumes a Salpeter-like IMF (ckimf = 2.35 which 
is consistent with that measured by Ferraro et al. 1997b). If the low-luminosity matter 
observed in the central 0.08 pc were entirely due to ss 1300 NSs of 1.4 M , then / ret ~ 30% 
(a reasonable value for collapsed clusters [Drukier 1996]). Alternatively, Ai/Cy ~ 6-7 can 
be explained by a Salpeter IMF if > 20% of the total population of heavy 1.0 M -WDs sank 
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into the NGC 6752 core during the cluster dynamical evolution. Either scenario must also 
be consistent with the observed shape of the star density profile. 

According to Cohn (1985), during the core collapse phase the surface density slope 
for the most massive component is expected to be a = —1.23 rather than the projected 
isothermal slope of a — —1.0. In the central cusp that forms during core collapse, the 
surface density slope of a component of stellar mass m is given approximately by: 

a = -1.89— + 0.65 

m d 

where rrid is the stellar mass of the dominant component. If the luminosity profile is domi- 
nated by turnoff stars of mass m = 0.8M Q and has a slope of a = —1.05, then the implied 
mass of the dominant, nonluminous component should be = O.89M , i.e. somewhat 
more massive than the adopted turnoff mass. This argument does suggest that the central 
gravitational potential is not likely to be dominated by a large number of neutron stars, but 
heavy white dwarfs still remain a possibility. 

Velocity dispersions provide a further constraint on the nature of the cluster potential 
well. The stars dominating the dynamics in the inner part of the cluster should have (see 
eq. 3.5 of Phinney 1992) a 1-dimensional central velocity dispersion a V:0 > 9-10 km s -1 . 
This is compatible both with the very wide published 2a interval for the a V:0 of NGC 6752 
(2.1-9.7 km s -1 , Dubath, Meylan & Mayor 1997) and with preliminary proper motion mea- 
surements of stars in the central part of NGC 6752, which would suggest a significantly 
higher one-dimension velocity dispersion a v>0 ~ 12kms _1 and the existence of strong ve- 
locity anisotropies (Drukier et al. 2003) 5 . Recent Fabry-Perot spectroscopy of single stars 
in NGC 6752 shows a flat velocity dispersion profile with typical one-dimension velocity 
dispersion a ~ 6-7 km s -1 within the central 1' (Xie et al. 2002). This is also marginally 
compatibile (given the ~ 10% uncertainty of the formula 3.5 of Phinney 1992). 

7.2. Case (ii): Local perturbator(s) 

We here explore the alternate possibility that the acceleration imparted to PSR-B and 
PSR-E are due to some local perturbator. 

NGC 6752 is a highly concentrated cluster and its core could host n* > 10 6 stars per 



5 If confirmed, such anisotropies would support the hypothesis of a relatively massive, probably binary 
black hole 
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cubic pc (assuming an average stellar mass of m ~ 0.5 M ). Hence close star-star en- 
counters are a viable possibility. In order to produce the line-of-sight acceleration |a;| = 
c\P/P\ = 2.9 x 10 _6 cms^ 2 seen in PSR-B (or PSR-E) a passing-by star of mass fn must 
approach the pulsar at s < (Gm\ai\~ 1 ) 1 / 2 = 0.0015 (m/0.5M Q ) 1 / 2 pc. An upper limit 
to the probability of occurrence of a suitable close encounter can be roughly estimated as 
3.7 x 10- 3 (m/0.5 M & f' 2 . Although this fi gure is not negliglible, the need of having 
two different canonical stars independently exerting their gravitational pull onto PSR-B and 
PSR-E, makes the joint probability of such configuration suspiciously low, < 10~ 5 . Rather 
than being randomly placed could the perturber be a binary companion to the MSPs? For 
such a companion not to have been aleady discovered by pulsar timing analysis would re- 
quire an orbital period Pb > 20 yr. Survival of such a wide binary is quite problematic in 
the core of a dense cluster like NGC 6752 — indeed no binary pulsar with Pb > 3.8 day has 
been detected in collapsed clusters to date. Seeing two such systems in NGC 6752 appears 
extremely unlikely. 

One may wonder if a single object, significantly more massive than a typical star in 
the cluster, could simultaneously produce the accelerations detected both in PSR-B and 
PSR-E. Recently, Colpi, Possenti & Gualandris (2002) suggested the presence of a binary 
black-hole (BH) of moderate mass (M bh+bh ~ 100-200 M ) in the center of NGC 6752 in 
order to explain the unprecedented position of PSR-A, a binary millisecond pulsar which is 
far away from the cluster center (D'Amico et al 2002). As shown in Figure 2, the projected 
separation of PSR-B and PSR-E is only d± = 0.03 pc. A binary BH of total mass M bh+bh , 
approximately located in front of them within a distance of the same order of d±, could be 
accelerating both the pulsars without leaving any observable signature on the photometric 
profile of the cluster. As the BH binarity ensures a large cross section to interaction with 
other stars, the recoil velocity v rec due a recent dynamical encounter could explain the offset 
position (with respect to C grav ) of the BH. However this scenario suffers of a probability 
at least as low as that of the previous two: placing the black hole at random within the 
core gives roughly a 1% chance that it would land in a location where it would produce the 
observed pulsar accelerations and the required t> rcc ^ 4-5kms _1 is also at the upper end of 
the expected distribution of v rcc for an intermediate mass black-hole in NGC 6752 (Colpi, 
Mapelli & Possenti 2003, in preparation). 



In summary, on the basis of the available data, the accelerations shown by PSR-B and 
PSR-E can be easilier accounted by the usually adopted hypothesis (i) (which could explain 
also the large positive P of PSR-D if it were not instrinsic), although the nature of the 
required extra amount (1500-2000 M & ) of low luminosity mass still remains puzzling. The 
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existence of local perturber(s) of the pulsar dynamics (case (it)) is a distinct possibility. 
While it seems extremely unlikely purely on the basis of the high value of P/P measured 
in PSR-B and PSR-E, there are other indications of a binary low-mass BH. In particular, it 
would explain: (a) the absence of any cusp in the radial density profile, (b) the flat velocity 
dispersion profile (Xie et al. 2002), (c) and (d) the ejection of PSR-A and PSR-C in the 
cluster outskirts (Colpi, Possenti & Gualandris 2002). However it is admittedly ad hoc 
hypothesis, requiring a fine tuned scenario. 

Clearly additional information must be collected. A longer baseline for timing measure- 
ments will allow to better constrain the presence of companions in very large orbits around 
PSR-B and PSR-E. In particular it will permit to derive (at least) upper limits on the second 
derivative of the pulsars' spin period, which is more influenced by by-passing stars rather 
than by the cluster potential well (Phinney 1992). Similarly, the (single massive or binary 
intermediate mass) black-hole hypotheses could be better investigated with spectroscopic 
determination of the dispersion velocity of the stars located in the pulsars' neighborhood, 
likewise with star density counts reaching higher magnitudes (hence exploiting a larger sam- 
ple of objects) and better spatial resolution (thus probing the innermost 1" of the cluster 
and the pulsars' surroundings). 
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Table 1. Surface brightness 
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rl and r2 are the inner and outer radii of each annulus in arcsec. 
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Fig. 1. — CMDs for stars in the two catalogs. Panel (a): The high resolution HST catalog in the 
(F555W, F336W - F555W) plane; Panel (b): The wide-field WFI catalog in the (V, B-V) plane. Only 
stars with r > 120" from the cluster center are plotted. Stars with V > 18.5 (dashed region) have been not 
considered in the density profile construction. 
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Fig. 2. — Computer map of the inner 20" x 20" region of the cluster. The large cross at coordinate (0,0) 
indicates the adopted C gra v The sma/( crosses are previous determinations of the cluster center: the C\ U m 
by Djorgovski (1993) is labeled with D and the C gra v by A089 (obtained considering only bright stars at 
V < 16) is labeled with AO. The nominal position of the 9 more central Chandra X-ray sources (Pooley et 
al. 2002) and of 3 MSPs are shown as empty star symbols and asterisks, respectively. The barycenter of the 
9 X-ray sources is also shown as a small cross and labeled with XRS. The dashed circonference encircles the 
core of the cluster. 
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Fig. 3. — Radial brightness profiles computed removing the stars brighter than V = 12, V = 13, V = 14, 
respectively, the faint limit always being V = 18.5 (see Figure 1). From the data of the lower plot, the best 
estimate of the central surface brightness is 15.07 mag arcsec~ 2 . Correction for extinction (using E(B — V) = 
0.04 , Harris 1996) leads to to 3.67 x 10 4 Ly© pc~ 2 . In the lower panel the data from A089 (open circles) 
are shown for comparison. 
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Fig. 4. — Observed radial density profile with respect to the the adopted C grav . The solid line is the best fit 
King model (r c = 13'.'7 and c = 2.1) to the observed density profile over the entire extension. Note that the 
adopted profile significantly undcrstimates the star counts in the region between 3" and 20". The adopted 
value for the parameter W (the central potential parameter defined by King 1996) is also reported. 
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Fig. 5. — Observed radial density profile with respect to the the adopted C gTav . The solid line is the best 
fit King model (r c = 28" and c = 1.9) to the outer points (that is, for distance larger than 10" from the 
cluster center). The King model has been combined with a constant star density of 0.25 stars arcmin~ 2 in 
order to account for the flattening of the density profile in the extreme outer region (r > 16.'5) of the cluster. 
The dashed line is the best fit King model to the 4 innermost points. 
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Fig. 6. — The central part (r < 1') of the radial density profile of NGC 6752. A modified power-law (from 
LCG95) has been used to fit the data. The solid line represents the best fit solution, whose parameters are 
also marked. The reduced %-squared is 0.28. 
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Fig. 7. — Maximum line-of-sight acceleration \ai max /c\ = \P/P\ versus radial offset with respect to the 
center of NGC 6752. The histogram represents the prediction based on the star density profile of §4 (nor- 
malized to the central surface brightness in the F-band obtained from the data of Figure 3) assuming a unity 
mass-to-light ratio. The dashed lines are analytical fits to the optical observations, labeled according to the 
adopted mass-to-light ratio. The measured values of P/P (filled pentagons, D'Amico et al. 2002) in the two 
MSPs with negative P (PSR-B and E) can be reproduced only for M/£v ~ 6 — 7. The open pentagons show 
our best guessed range of maximum \ai/c\ for PSR-D: the upper value is calculated assuming a negligible 
intrisic positive P sc i', the lower value is estimated taking into account intrinsic P sc i from the observed scalings 
between X-ray luminosity and spindown power for MSPs (see D'Amico et al. 2002 and reference therein). 
Given the relative large uncertainty, the value of \ai/c\ it is not further taken into account in the discussion. 



